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1. Introduction 

The formulation and application of the string hypothesis has been a story of success as well 
as a longstanding mystery in the studies of exactly solvable models [1,2,3,4]. For finite sys- 
tems, the deviations from this hypothesis have been discussed [5,6,7,8] and prove essential 
for a proper description of critical properties [9,10,11]. On the other hand, for finite tem- 
peratures we expect the hypothesis to give the correct free energy in the thermodynamic 
limit [3,4]. 

An alternative approach to thermodynamics which is free of the string hypothesis has 
been proposed [14,15,16,17,18,19,20,21,22,23]. There the essential idea is to deal with 2D 
classical counterparts. A novel object, "the quantum transfer matrix" plays a fundamental 
role. Remarkably, the problem of summing up all states reduces to the evaluation of 
only the largest eigenvalue of the quantum transfer matrix. Of course, the remaining 
task is the diagonalization of the quantum transfer matrix. In order to achieve this, a 
more sophisticated approach has been proposed recently by incorporating the classical 
integrability structure in two dimensions [20,21,24,25,26,27]. The difficulty still remains in 
the fact that the resultant system is virtually of finite size with interactions depending on 
the size. 

Locating the Bethe ansatz roots for finite system sizes is another difficult task, which 
we avoid here by one of the following two alternative strategies. We encode the information 
of BAE roots for finite size systems by 

(1) a finite number of coupled non- linear integral equations for a finite number of unknown 
functions derived within a direct Bethe ansatz for the diagonalization problem, or 

(2) an infinite number of coupled non-linear integral equations for an infinite number of 
unknown functions obtained within the fusion approach to the diagonalization prob- 
lem. 

Note that both transformations are exact. However, approach (1) has an evident advan- 
tage in numerical investigations. Actually several thermodynamical quantities have been 
explicitly obtained for extensive ranges of temperature and chemical potentials for ID cor- 
related electron systems within this scheme [25,26,27]. The latter approach (2) would be 
less important in actual calculations. Nevertheless, it is interesting in its connection with 
the traditional string approach. The infinite number of unknown functions introduced 
there, physically corresponds to the fused ( "higher spin" ) transfer matrices. These trans- 
fer matrices satisfy algebraic functional relations. For the simplest case, these relations 
are shown to have a deep mathematical origin, an exact sequence of Yangian modules 
[28,29,30,31]. 

In this paper we will show that these functional relations can be transform,ed into non- 
linear integral equations, which are identical to TBA equations. Note that our approach is 
completely independent of the string hypothesis. It thus gives an independent derivation 
and further support for the validity of the resultant TBA equations. For the physical 
applications we need only the eigenvalue of the fundamental ( "un- fused" ) quantum transfer 
matrix. Thus the introduction of the whole commuting family of transfer matrices seems 
to be redundant. The point is that the exact relation among them makes the evaluation 
of the fundamental one possible. This fact was first realized in [20] for the RSOS chains. 
In this paper, we will apply this scheme to two ID highly correlated electron models to 
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demonstrate that the relation of the fusion analysis of the quantum transfer matrix and 
the TBA is not accidental but universal. 

The novel formalism, utilizing the functional relations, not only reproduces the results 
obtained from the string hypothesis, but goes well beyond. TBA equations usually deal 
with the equilibrium free energy, i.e. we can only deal with the "largest eigenvalue sector" 
in the sense of the quantum transfer matrix. The functional relations are, on the other 
hand, valid for arbitrary eigenvalue sectors. The latter, therefore offers a wider possibility 
in calculating various physical quantities. Actually, we need to evaluate the subleading 
eigenvalues for obtaining correlation lengths. Recently, some attempts have been made in 
the context of deformed conformal field theories, to modify TBA equations so as to deal 
with excited states [33,34,32,20,36,37,38]. We will show how the quantum transfer matrix 
formalism naturally leads to excited state TBA equations. To our knowledge, this is the 
first explicit derivation of the excited state TBA equations for ID lattice electron systems. 

This paper is organized as follows. In the next section, we will give a brief review of 
the quantum transfer matrix approach to one- dimensional quantum systems at finite tem- 
peratures. In section 3, the functional relations among the relevant transfer matrices are 
given. The corresponding relations for the ordinary row-to-row transfer matrices were ob- 
tained recently [39] . Here, we will summarize the necessary modifications for the quantum 
transfer matrices. Section 4 and 5 are devoted to the thermodynamics of the integrable 
t — J model and the supersymmetric extended Hubbard model, respectively. We will show 
explicitly how the functional relations of section 3 can be transformed in the case of the 
largest eigenvalue into the TBA equations of [40,41]. To this end, we adopt assumptions 
about the analyticity of auxiliary function, which are checked numerically for the case of 
finite Trotter numbers. In section 6, we will discuss "excitations at finite temperatures", 
i.e. the correlation lengths of static two point functions. Assisted by numerical investiga- 
tions, the excited state TBA equations are easily derived there. Section 7 is devoted to 
the summary and discussion. 
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2. The quantum transfer matrix formalism 

In the foUowing we wiU mainly focus on two specific models of ID highly correlated electron 
systems: the supersymmetric t — J model [43,40,44,45] and the supersymmetric extended 
Hubbard (SEH) model [42]. They describe spin-1/2 electron systems on the lattice. The 
prototype of such systems is the celebrated Hubbard model. The supersymmetric t — J 
model may be viewed as the large Coulomb repulsion limit of a generalized Hubbard model. 
The SEH model on the other hand, generalizes the Hubbard model by additional bond- 
charge interactions etc. It shares the same interesting physical symmetry with the standard 
Hubbard model resulting in the eta pairing [46,47]. Mathematically the supersymmetric 
t — J model and SEH model are in the category of models based on Lie superalgebras gl{r\s) 
with r = 2, s = 1 for the t — J model and r = s = 2 for the SEH model. As remarked 
in [50,51,24], their classical counterparts are special cases of the Perk-Schulz [48,49] model 
with rational vertex weights 7?.(j^^ 

K^(t;) = l + e«^, 

nZM^e^e^v, (2.1) 

which satisfy the Yang-Baxter equation. As a consequence the row-to- row transfer matrix 
T{v) constitutes a commuting family with respect to the spectral variable v. Cq = ±1 
are discrete parameters determining the 'grading' of the system. The Hamiltonian of the 
associated quantum system is given by the logarithmic derivative 

7Y = -^lnr(^)|,=o. (2.2) 
dv 

The supersymmetric t ~ J Hamiltonian is obtained for (ei,e2,e3) = (-|-,— ,-|-). Likewise 
the grading (ei, 62, 63, 64) = (+, — , — , +) leads to the SEH model. 

Next we introduce Boltzmann weights and related to by anti-clockwise and 
clockwise 90^ rotations 



nZiv) = nf,{v), n^iv) = K^{-v). (2.3) 

According to (2.1), (2. 2) we find 

T{v) = Tfle^^+^('''), T{v) = ri,e-^+<^('^'), (2.4) 



where 7r ^ are right and left-shift operators. Therefore the partition function of the 
quantum system is given by [17,20,22] 

Z^lim Ty[T{u)T(u)]^^\ u^-(3/N. (2.5) 

The right-hand side of this equation may be viewed as the partition function of a staggered 
Perk-Schulz model consisting of alternating rows of TZ{u) and TZ{u) operators. For the 
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following analytic calculations the column-to-column transfer matrix T^'^^ (quantum 
transfer matrix) is best adapted. For quite general systems it can be shown that T^'^^ 
possesses a gap between the largest and next-largest eigenvalues A^ax and A^ext persisting 
in the limit — > cxd. Therefore, the problem of computing the free energy / per site of 
the quantum system at finite temperature is reduced to just the evaluation of the largest 
eigenvalue Amax 

/ = -/cBTlnA„,ax, (A^^oo). (2.6) 

If also the gap can be calculated (as shown below) the correlation length at finite temper- 
ature can be derived 

e=(ln^) \ (AT^oo). (2.7) 

The quantum transfer matrix is given by an alternating product of 7l{u) and TZ{—u) 
operators. The important observation [20,52] is that TZ{v) and 7^(f ) share the same inter- 
twiner as can be proven most easily in a graphical way, see Fig.l. Therefore 

N/2 ® 

T^™{u,v)= W n{v + u)®n{v-u), (2.8) 

represents a family of commuting matrices comprising the 'physical' QTM at f = 0. This 
intertwining property implies that the Yangian algebra, in the realization RLL = LLR, is 
identical for both monodromy matrices. Consequently, we should have the same functional 
relations among the transfer matrices in spite of apparent differences between the explicit 
eigenvalues. Of course this argument needs further elaboration. However, we leave it as an 
interesting future problem and assume the validity of the above argument in the following. 
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3. The fusion hierarchy- 
First, we like to sketch the strategy for the diagonahzation of the quantum transfer matrix 
for the Perk-Schulz model spectral parameters u, iv (where we have introduced the factor 
i for later convenience) . In the following analysis we will be dealing with transfer matrices 

formulated for fused Boltzmann weights. These matrices arc denoted by Tjn\u,v) where 

Of course T^^\u,0) is the only quantity we are interested in. 

For calculating the eigenvalues of T^^\u,0) it is essential to deal with all Tm\u,v) as 
our method employs functional equations dealing with the dependence on v and all other 
matrices for any a, m [20] . 

In order to present our results for the eigenvalue of the quantum transfer matrix we 
adopt a compact notation using the "Yangian analogue of Young tableaux" [53,54,55]. 
First consider the simplest case: the quantum transfer matrix with r + s dimensional 
degrees of freedom on both vertical and horizontal edges. Let us introduce boxes with a 
letter e {1, • • • , r + s}. Each box corresponds to an expression as follows. 

^y„(....K(.-,.)e---^-- '^°^--°"^°^;-°' , (3.1) 



where 



fa{x) 
9a{x) 



(1 + €ix)^/^, for a = 1, 
x^/'^, otherwise, 

(1 + €r+sx)^/^, for a = r + s, 
x^/^, otherwise. 



Q^(x) = I n7=i(^ - xf), for a = 1, ■ ■ • , r + s - 1, 
[1, for a = 0, r + s. 

Note that a box carries the spectral parameter dependence. In the last equation, x^"'^ {j = 
l,---,Na) denotes the Bethe ansatz roots with "color" a. (We adopt the convention 
Nq = Nr+s = 0). As derived by an algebraic Bethe ansatz the eigenvalue A^\u,v) of the 
quantum transfer matrix is of the form 

1 — \-s 

A^^\u,v) = Y,[^. (3.2) 



The Bethe ansatz equations are obtained by the pole free condition for the eigenvalue. 
The BAEs are therefore given by 



Res (a) ( a + a + 1 ) = 0, 

v=x) It, I — ! Iv^ ' 



ta+l[u + ix] ')ga+l[u- IX) ') 

Qa-i{.xf^)Qa{xf^ - iea+i)(5a+i(a;J"^ +iea+i) 

Qa-l{xf^ - iea)Qa{x'f^ + iea)Qa+l{x^^'') 



(3.3) 
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Starting from A^^''^'' as in [54] , one can generate a set of analytic functions under the BAE. 
The resultant functions are expected to be eigenvalues of fusion transfer matrices. For the 
simplest case like s/2 spin chains, we can explicitly prove this [56,57]. We do not attempt 
to prove this observation for the general case as our analytic treatment below does not 
depend on it. In the case of row-to- row transfer matrices, the "Bethe-strap procedure" 
has been executed for sl{r\s) models [39]. The eigenvalues of fusion transfer matrices 
are parameterized by a similar set of tableaux under certain combinatorial rules. A set 
of functional relations among them is also found. As remarked above, the same set of 
functional relations holds for the quantum transfer matrices as well as the combinatorial 
rules for the tableairx giving the analytic eigenvalues for the fusion transfer matrices. If one 
replaces boxes for the row-to-row case by the quantum ones in eq (3.1), then the explicit 
eigenvalues of "fused transfer matrices" are obtained. Let us further explain the rules in 
the present context. Consider an a x m rectangular tableaux. We fix the "coordinate" so 
that the upper-leftmost box is in the position (1,1), the one box to the right is in (1,2) and 
the one box lower is in (2,1). We associate the spectral parameter v + i{m — a) /2 + i{j — k) 
with the box in {j, k). Denote the letter in the (j, k) box by £j^k- We identify such tableaux 
with the expression 



nn 

j=l k=l 



v-\-i{m — a) /2+i{j — k) 



(3.4) 



To describe the combinatorial rules we introduce the order among letters, 1 -< 2 -< • • • ^ 
r+s. These letters are further classified into two subsets J+ and J_. Then the admissibility 
conditions are [39] 

(1) Ij-i^k -< ij,k if £j,k e J+ 



(2) 
(3) 
(4) 



^j — l,k — ^j,k 
^j,k — l — ^j,k 
1 ^ " 



Define A^^ 



if ij^k e J_ 
if £j^k e J+ 
'j,k if £j,k e J- 

{u, v) by 



{ij^k} under rules (l)-(4) J=ifc=i 



^j,k 



v+i{m—a)/2+i{j — k) 



(3.5) 



(-u, v) is expected to be the eigenvalue of the quantum transfer matrix whose trace has 
been taken over the Yangian module Wm'' ■ We have one further rule [58] : 

AI^\u, v) = if a > r and m > s. (3.6) 



The functional relations read 

^tH^,v-'-)Al^Hu,v+'-)^Ai:+'Hu,v)Ai^-'\u,v)+Al^U^^ (3.7) 

with the convention Aq"'^ = A^^ = 1. These relations can be proved using a quantum ana- 
logue of the Jacobi-Trudi formula and Plucker's relation [55,39]. Recently some attempts 
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have been made in the study of finite-size corrections int — J like models employing similar 
functional relations [59] . There the row-to- row transfer matrix was considered and a subset 
of fusion functions was used. 

In the next two sections, we will transform the complete set of functional equations 
into integral form which is useful for further analytic treatment. We will show explicitly 
how the simple functional relations above are nothing but TBA equations for the quantum 
chains related to the Perk-Schulz model. We will focus on two specific examples of interests, 
the integrable t — J model and the SEH model. 

4. Fusion based derivation of TBA for the supersymmetric t — J model 

For the t — J model, we choose a convenient grading ei = €3 = l,e2 = —1. Then J+ = 
{1,3}, J_ = {2}. For example, we have the following five tables for A^^^ 



11, 12, 13, 23, 33, 



and four tables for A^-*, 



1 




1 




2 




2 


2 




3 




2 




3 



and so on. From the analytic point of view, it is better to divide the A^"* (w, v) as defined in 
(3.5) by a common factor, such that the resultant expression is a polynomial of degree 
in V irrespective of the values of a and m. For that purpose we prepare the normalization 
functions, 



m — 1 



(m + 1) . 



{ni + 1) , ., 



a 



n'^^\x)=\{h'i\x+{3-^^)i)x 



i=i 

a-l 

n 



(4.1.a) 



(4.1.6) 



where 4>±{v) = (v ± iu)^/'^. Then the normalized functions A^^ are given by 

A(,«)(«,x) = n(;j)(a;)AL"^(a;), 



h.^^{x) = (f)+(x- ^i)(t)-(x + yz), 
A(«)(x) = 0+(x+|i)0_(x-|i). 



(4.2.a) 
(4.2.6) 

(4.2.C) 



Hereafter we drop u from Am"* for simplicity. 
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Now the functional relations read, 



A^rk\x - V2)A^^(a; + i/2) = A^^i(a;)A^Li(a;) + A^^(a;)A^)(a;), m = 1, • • • , oo,(4.3.a) 
A^^\x-i/2)A^^\x + i/2)=A^^\x)Al^\x)+A^^\x)A^f\x), (4.3.6) 
A^^\x-i/2)A^^\x + i/2)=A^^l^{x)A^^l^{x), m = 2,---,oo, (4.3.c) 
A[''\x-i/2)A[''\x + i/2) = A[''+^\x)A[''-^\x), a = 3, •••,00. (4.3.d) 

In the above, not all Am\x) are independent; we have the identity A^^'^^\x) = A^\x), 
which can be easily proved by explicit forms. Now we apply transformations generalizing 
the ones in [20,57,31] 

Yi'Hx) = ^^^^M^^, form>l, (4.4.a) 

A^^(a;)A^)(a;) 

= for a > 2. (4.4.6). 

A["-^\x) 

We call the functional relations in terms of the functions (x) simply the "Y-system" 
[12,35,60]. The explicit Y-system is now given by 



{1 + {Yi'\x))-^y 

Yi'\x - i)Yi'\x+'-) = {1 + Yi%(x))il + Yi%ix)), for m > 2, (4.5.6) 

y(2) , i . (2) . i Mx + V2)(^- {x - i/2) y/'^ jx) 

^^+2^-0_(x + z/2)</>+(x-z/2)(i + ya)(^))' ^^"^-^^ 

Yi'\x '-)yI'\x + ^) = , (4.5.d) 

' ^ 2^ ' ^ 2^ (l + (Ff)(x))-i)' ^ ^ 

yI"\x - '-)yI''\x + |) = Y'/"+')(x)r/"-')(x), for a > 4. (4.5.e) 

At this stage we note some numerical results concerning the analytic properties of Am\x). 
They are analytic due to the BAE, and have AT zeros. Through finite AT studies, we^find 

their characteristic behavior. For the eigenstate giving the largest eigenvalue for A\^\x), 
the zeros are located on curved lines with imaginary parts close to ±(m+l)/2. Therefore all 
Ym\x) except for a = m = 1 are Analytic Nonzero and have Constant asymptotics(ANZC) 
in the strip ^x e {—1/2, 1/2}. Similarly we have numerically checked that 1 + Ym^ (x) has 
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the ANZC property in this strip. For the case of (1 + 1^2^^^) see Fig. 2. To cancel zeros and 
holes for y}^\x), we define 



{ {tanh^{x + i{l + u))t8inh^{x-i{^+u)))-^/^Yjh\x), for a=m=l, 

\^m"'*(a;), otherwise. 

_ (4.6) 

Now the left hand sides of eqs. (4.5) are invariant if we replace Ym\x) by Ym\x). After 
this replacement, both sides of eqs. (4.5) are ANZC functions in a narrow strip including 
the real axis. Therefore the Fourier transformation of logarithms of both sides can be 
simplified due to Cauchy's theorem. This results in a set of linear equations in Fourier 
space which can be solved easily. The reverse Fourier transform leads to the following 
infinitely many coupled non-linear integral equations. 

logYl'\x) = ei{u, x)+Go* log(l + Y^^'^){x) - Go * log(l + ^)(^), (4.7.a) 

logYi'\x) = Go * log(l + Y^^l,){x) + Go * log(l + Y^^l,){x), for m > 2, (4.7.6) 
logy/^^(a;) = e2iu, x) + Go * logYl^\x) - Go * log(l + yI^^){x), (4.7.c) 
logyf ^(x) = Go * logYl^\x) - Go * log(l + ^), (4.7.d) 

log y}"^ (x) = Go * log y/""^^ (x) + Go* log Y"/"+^^ (x), for a > 4. (4.7.e) 
Here A* B{x) means convolution A* B = J A{x — y)B{y)dy, and 

/ TT 1 TT 1 \^/2 

ei{u, x) — log! tanh —{x + i{- + u)) tanh —{x — i{- + u)) j , 

r e~l'=l/2-*'=^^ sinh/cM 
2fccosh/c/2 
1 r e-*'^^ 1 



27r 7 2coshA;/2 2cosh7ra;' 
Note that logy/^^(a;) can be solved in term of log(l H — from the last two relations in 



the above. Resubstitution of the result into the third relation above leads to 

log Y-/'^(x) = e2{u,x) - Gi * log(l + -^){x) - Go * log(l + Y'/'^)(x), (4.7.cO 

where Gi{x) = ^ j dk^^^^^. 

Therefore we have the resultant Y-system including Ym \ (m = 1, • • • , oo) and Y^'^\ The 
quantity of our interest, the free energy, is given in terms of Ym^ as, 

-/?/= hm logAi^)(« = -^,0)= hm log0+(i)0_(-i) - log Y-/')(« = 0) (4.8) 

AT— >oo iV N—K30 I\ 
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Now the limit A?" — >^ oo is taken easily by analytic means. To write down the final result, 
we make slight changes of notations 



Yl'\x)=exp{-e{x)/T), 
Yi'\x)=exp{<j)m-i{x)/T), 
yI^\x) ^exp{-ij{x)/T). 



m > 2, 



Then we have 



e = 277^0 + TGq * log 



(l + e-^i/^)' 

iP = 27rGi + TGo * log(l + e"^/^) + TGi * log(l + e'^/'^), 
01 = TGo * log[(l + e-^/^)(l + e-^^/^)], 

= TGo * log[(l + e<^--/^^')(l + 6-^-+^/^)], 
/ = 1-V(0). 



(4.9.a) 
(4.9.6) 
(4.9.C) 



(4.10.a) 

(4.10.6) 
(4.10.C) 
(4.10.d) 



They are nothing but the TBA equations obtained in [40] for Xs = 1,A' = 0. The shift 
by unity in the free energy is due to a trivial shift in the definition of the Hamiltonian. 
Thus we have succeeded in deriving these equations independent of the string hypothesis. 
In the next section, we further apply the same strategy to another model in ID correlated 
electron systems to convince ourselves of the universality of the result. 

5. TBA for the SEH model 

For this model, we choose the grading {ei, £2, €3, £4} = {+,—,—,+}. Then J_|_ = {1,4}, 

J_ = {2, 3}. Due to these choices, we have non- vanishing A's: A^rk\ A^\ A^^^ A^"\ (m = 
0, • • • , 00, a = 0, • • • , 00). While the derivations of TBA based on the string hypothesis for 
the t — J and the SEH seem to be different, we find them quite similar by adopting the 
present strategy. We thus sketch the relevant definitions and relations for the SEH model, 
without going into details. As for the t — J model, we re-normalize Am\u, v) by the same 
nm\iJ',v) in (4.1.b). The renormalized functional relations are given by 

' {x + i/2) = A« , (x) A (x) + a(^) {x)A^^^ (x) , 





(x 


-i/2)A^) 




{x 


- V2)Ai") 


A^ 


{x 


- i/m? 


A^ 


{x 


- V2)A^^ 


A (a) 
^^2 


{x 


- V2)A^"^ 



m 



{x + i/2) ^ Ai''\x)A\;;'{x)+A\''~'\x)A\^^'\x), 



(a) 



a = 



(x + i/2) = Af(a;)ArV) + A'^'\x)A'^' (x), 

(x + z/2)=A^^,(x)a2^(x), 



(a; + V2) = A^"+'^ 



{x)Ai^-'\xl 



m — 3,- 
a = 3, 

In addition we have the equalities, 

A^^\x) = A'^^\x), m>2. 
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00, 



, 00. 



■ , 00, (5.1. a) 
2,---(6ol,.6) 
(5.1.C) 
(5.1.d) 
(5.1.e) 

(5.2) 



We introduce the Y functions as follows 



y}''\x) 



A^\x)A^{x) 

A^r^\x)Ag\x) 
^\x)^Hx) 
A^^^^\x)A^r^\x) 

Then the following Y-system holds, 

{l + Yi'\x)) 



(1) 



for m > 1, 



for m = 1, 2, 



for a > 3. 



(5.3.a) 
(5.3.6) 
(5.3.C) 



(l + (Ff)(a;))-i)' 

(l + yi^)(a;))(l + y/^)(a;)) 

(i + (yi')(x))-i) ' 

{1 + yP{x)) 



for m > 3, 



(i + (yf)(x))-i)(i + (yr(^))-^) 



(1)/ 



y/"+')(x)yi^"-^^(x) 



(a-l), 



for o > 3, 

(2), 



(t)+{x + z)(/)-(a; - i) (1 + (yi^'^(x))-i) 
<t)^{x + i)<j)+{x - i) {i + Y^^\x)) 



(5.4.a) 

(5.4.6) 
(5.4.C) 
ibA.d) 
(5.4.e) 
(5.4./) 



Numerically, both sides of the above equations have the ANZC property in a narrow strip 
near the real axis, except for the Ihs of the first equation. As for the t — J model, we 
cancel zeros and singularities by multiplying the scalar as in (4.6). Then the Y-system is 
equivalent to the following set of non-linear integral equations 



logy/'^(a;) = ei(«,x) + Go * log(l + Y^'>){x) - Go * log(l + ^){x), 

Yi 



1 



(5.5.a) 
(5.5.6) 



(i + (yf)(x))-i) 

logYi'Hx) = Go * log(l + Yi'l,){x) + Go * (1 + YX';,){x)), for m > 3, (5.5.c) 



log y/"^ (a;) = Go* logy/" ^^(x) + Go * logy/"+^^(a;), for a > 3, {5.5.d) 
logyi')(a;) = e3{u,x) + Go * log(l + {yI''\x))-^) - Go * log(l + Yi^^){x), (5.5.e) 
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logFf = Go * log[ .3/^ + ^"^''^^"^^ ], 

{1 + {y}'\x))-^){1 + {y}'\x))-^) 

= e2{u,x) + Go * log[ (l + (y2 (^))"') ], (5.5./) 

{i + {yI'\x))-^){i + yI'\x)) 

e3{u,x) = Nj 2/,cosh/c/2 
A^"*^^ (tt, a;) in terms of yiT^ reads, 

00 

A\'\u,x) = log(/)+(x + (x + log(l + (y/'y (5.6) 

Now let — > cxD under identifications, 

a„ = yW, (n>l); K = 1), 

Is = l/Y^'^. 

Then the resultant TBA equations are nothing but the ones in [41], after taking some 
appropriate linear combinations 

log CKn = +^0 * log[(l + an+l)(l + "n-l)] 

- 5^,i(27rGo/T + log(l + l//3i) - 5^,2 log(l + l//?2), n > 1, 
log 7, = +Go * log[(l + 7s+i)(l + 7,_i)] + d,,iGo * log(l + l/A), s > 1, 

(l + /5i) 1 



log A = -27rG'o/T + Go*log[ 



(l + cti)(l+7i)- 



log/?2 = -27rG'i/T + Go * log[ ^) + ^^^;^ ], (5.7) 

(1 + 0:2) 

f=l-J2Kn*log{l + l/an). 
n 

Note that the present choice of Hamiltonian corresponds to fi + U/2 = 0. 
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6. Excited state TBA 

So far we have shown that the quantum transfer matrix formahsm for the largest eigenvalue 
of A^^^"* with adequate information on analyticity properties directly relates the functional 
relations with the TBA equations. The investigation of other eigenvalues naturally leads 
to the excited state TBA equations, which characterize correlation lengths at finite tem- 
peratures. Note that the functional relations are valid irrespective of the sectors under 
consideration. The only difference to the "ground state" is the modification of the analyt- 
icity conditions by additional zeros or singularities entering into the strips where auxiliary 
functions in the ground state are strictly analytic and non-zero [57,61]. Let us examine 
explicitly the spin excitation for the solvable t — J model. 

For the ground state, we have ^x''^ and ^x'^^'' Bethe ansatz roots complex conjugate 
to each other. This time, we also have such complex conjugate pairs of BAE solutions, 
however, the number of solution is N/2 — 1 for each. This does not seem to be a serious 
modification of the ground state, however it results into considerable differences in the zeros 
for A^'*- We find N — 2 zeros for A^"* on slightly curved lines close to "^x ~ ±(m + l)/2 
as for the ground state. Note that two zeros are missing on these lines. Instead we find 
(by numerical investigations for vanishing chemical potential) the additional structure: 

(1) A^-* has zeros at x^^_^ on the real axis. The strictly largest eigenvalue is given by two 
symmetrically distributed zeros x^j_ = ±Xm^ . 

(2) A[^\ a > 2, has zeros at ±(0;^'' + {a — 2)i/2) on the imaginary ax;is. We find x^^^ is 

(2) 

pure imaginary and \x\ \ > 1/2. 

As a corollary from the identity A^'' = A^™'"'"^'* we get 

(3) Am has zeros at ±{xi ' + {a — l)i/2) on the imaginary axis. 

Note that the behavior of zeros (2) is consistent with the "trivial" T-system, 

a[''\x - i/2)A^f{x + i/2) = A^^+^\x)A^^-^\x), a = 3, • • • , 00. 

Thus, in terms of the Y functions, we find the following zeros and singularities in the strip 
'^x e {-1/2,1/2}, 

Yi^\x) = 00, at X = X^^^. (6.1) 

In order to have ANZC functions in the strip we must modify the Ym^ functions of section 4 
further to 

yS^^Y^^\{ tanh |(x - J tanh |(a: - x^.,, J, 

cr=± 

lf)^lfVntanh|(x-x«). (6.2) 
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The resultant are now ANZC functions in the strip. Repeating the same argument 
as in section 4, we thus arrive at 

logYl^'{u,x) = ei(u,x) + ^0 *log- ^ (x) + ^ log(tanh -(x - 4 a))' (6.3.a) 

(1 + 1/^1 ) a=± ^ 

logy a;) = Go * log(l + Y^^\){\ + y'^1^)){x) 

+ X] log(tanh -(x - xj^^i,^) tanh - x^^_^^^)), for m > 2, (6.3.6) 

logyf = e2(«,x) - Gi *log(l + ^)(x) - Go *log(l + 

- log(tanh ^(x - x^l\)), (6.3.c) 

and the x^m.a satisfy, 

rWKx = a;^,,±V2) = -l. (6.4) 

In the new integral equations the convolutions are still evaluated with integral contours 
along the real axis. In this way, we have derived quite naturally excited state TBA equa- 
tions for spin excitations. We want to stress once more the meaning of these excitations. 
The corresponding eigenvalue determines via (2.7) the decay of static correlation functions, 
it is not related to dynamical correlations or "excited energies at finite temperature" [1]. 
For other types of excitations, one needs to identify the patterns and locations of zeros and 
singularities of the auxiliary functions. This is possible for small Trotter numbers, as we 
know explicitly the spectrum from "brute force" numerical calculations. We expect the so 
obtained patterns would be valid also in the limit N ^ oo. Lastly, we want to comment 
on our derivation of the TBA equations on the basis of ANZC properties. It appears that 
any solution to the integral equations is a solution of the functional equations and thereby 
yielding an eigenvalue to the quantum transfer matrix. The question remains whether the 
obtained eigenvalues are the largest ones. This however, is a byproduct of this section. 
Any time we modify the ANZC properties by allowing zeros to enter the physical strips the 
eigenvalue is lowered. This is the analytic justification of the previously employed analytic 
properties of the largest eigenvalue. 
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5. Summary and Discussion 

In this paper, we have shown the equivalence of the fusion relations among the quantum 
transfer matrices and the TBA equations under appropriate assumptions about analytici- 
ties of auxiliary functions. This has been demonstrated by adopting two examples of ID 
highly correlated electron systems. 

In our approach the derivation of the TBA equations is based essentially on ANZC 
properties of the auxiliary functions satisfying the functional relations. In some sense, 
the string hypothesis is now replaced by ANZC assumptions. We like, however, to point 
out the following. First of all, the analyticity properties are easier to control. We have 
checked numerically the ANZC properties for various configurations and for several Trotter 
numbers up to A?" = 32. Second, in contrast to deviations from the string hypothesis for 
the finite system, we do not observe numerically any violation of the latter assumption for 
finite Trotter numbers. 

Another advantage in our approach is the straightforward derivation of "excited state 
TBA". The additional chemical potential terms in TBA equations in the CFT limit of 
previous studies [33,34,35,32] is now replaced by the concrete objects, the contribution by 
the additional zeros Xm^ of the fusion transfer matrices. 

For the actual calculation of physical quantities such as the specific heat or correlation 
lengths, we must deal with, in principle, an infinite number of unknown functions and 
parameters x^'^ . Let us sketch the situation in the case of the t — J model. In order 
to solve numerically the TBA equations (4. 7. a)- (4.7.e) we truncate the infinite set of 
integral equations (4.7.b) to a finite one by approximating logYm\x), for m ^ 1 by its 
asymptotic value. We thus introduce the quantity log{Y^\o) / Ym\(X))) as an estimate 
for the systematic error. To evaluate this, we use the explicit solutions to BAE for finite 
Trotter numbers N together with functional relations, thereby allowing the evaluation of 
these quantities with high numerical precision. In the limit m ^ 1 we find the following 
asymptotic behavior 

log(FW(0)/yW(oo))~C/m^ 

For the example = 12, P = 1/T = 3.5 (and vanishing external fields), we have C ~ 6. 
Thus we need more than 2500 auxiliary functions Ym^ to obtain a numerical error less 
than IQ-^. This already may be a sign of difficulties in obtaining good numerical accuracy 
for the traditional TBA analysis. 

Fortunately, as we already remarked in the introduction, we have an alternative ap- 
proach where only a finite number of auxiliary functions is necessary [21,24,25,26,27]. Al- 
though the connection to the traditional TBA is obscure, this formulation is more efficient 
in numerical investigations for physical quantities at finite temperatures. For the t — J 
model for instance, we need only three auxiliary functions for a closed set of non-linear 
integral equations. This results in a numerical accuracy better than 10~^ [25,26] for the 
specific heat, compressibility and other thermodynamic quantities. Now the evaluation of 
correlation lengths is under investigation. We hope to report on explicit results based on 
this formulation in the near future. 
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